Testing satellite and ground thermal imaging of low-temperature fumarolic fields: The dormant Nisyros Volcano (Greece)

The Nisyros Volcano (Greece) was monitored by satellite and ground thermal imaging during the period 2000–2002. Three night-scheduled Landsat-7 ETM⁺ thermal (band 6) images of Nisyros Island were processed to obtain land surface temperature. Ground temperature data were also collected during one of the satellite overpasses. Processed results involving orthorectification and 3-D atmospheric correction clearly show the existence of a thermal anomaly inside the Nisyros Caldera. This anomaly is associated mainly with the largest hydrothermal craters and has land surface temperatures 5–10 °C warmer than its surroundings. The ground temperature generally increased by about 4 °C inside the main crater over the period 2000–2002. Ground thermal images of the hydrothermal Stephanos Crater were also collected in 2002 using a portable thermal infrared camera. These images were calibrated to ground temperature data and orthorectified. A difference of about 0–2 °C was observed between the ground thermal images and the ground temperature data. The overall study demonstrates that satellite remote sensing of low-temperature fumarolic fields within calderas can provide a reliable long-term monitoring tool of dormant volcanoes that have the potential to reactivate. Similarly, a portable thermo-imager can easily be deployed for real-time monitoring using telemetric data transfer. The operational costs for both systems are relatively low for an early warning system.     

Universal time from VLBI single-baseline observations during CONT08

The IVS Intensive sessions are single-baseline, 1-h VLBI sessions carried out everyday in order to determine Universal Time (UT1). We investigate different possibilities to improve the results of such sessions. We do this investigation by extracting 2-h single-baseline sessions from the CONT08 data set. These are analysed like normal Intensives, and the results are compared to the results of the analysis of the full CONT08 data set. We find that tropospheric asymmetry is the major error source for the single-baseline sessions. It is possible to improve the accuracy of the estimated UT1 either by using accurate a priori tropospheric gradients or by estimating gradients in the data analysis.     

Measuring tree stem diameters using intensity profiles from ground-based scanning lidar from a fixed viewpoint

This paper presents a method for using the intensity of returns from a scanning light detection and ranging (lidar) system from a single viewing point to identify the location and measure the diameter of tree stems within a forest. Such instruments are being used for rapid forest inventory and to provide consistent supporting information for airborne lidars. The intensity transect across a tree stem is found to be consistent with a simple model parameterised by the range and diameter of the trunk. The stem diameter is calculated by fitting the model to transect data. The angular span of the stem can also be estimated by using a simple threshold where intensity values are tested against the expected intensity for a stem of given diameter. This is useful when data are insufficient to fit the model or the stem is partially obscured. The process of identifying tree positions and trunk diameters is fully automated and is shown to be successful in identifying a high proportion of trees, including some that are partially obscured from view. The range and bearing to trees are in excellent agreement with field data. Trunk angular span and diameter estimations are well correlated with field measurements at the plot scale. The accuracy of diameter estimation is found to decrease with range from the scanning position and is also reduced for stems subtending small angles (less than twice the scanning resolution) to the instrument. A method for adjusting survey results to compensate for trees missed due to obscuration from the scanning point and the use of angle count methods is found to improve basal area estimates and achieve agreement within 4% of field measurements.     

VLBI-derived troposphere parameters during CONT08

Time-series of zenith wet and total troposphere delays as well as north and east gradients are compared, and zenith total delays (ZTD) are combined on the level of parameter estimates. Input data sets are provided by ten Analysis Centers (ACs) of the International VLBI Service for Geodesy and Astrometry (IVS) for the CONT08 campaign (12?C26 August 2008). The inconsistent usage of meteorological data and models, such as mapping functions, causes systematics among the ACs, and differing parameterizations and constraints add noise to the troposphere parameter estimates. The empirical standard deviation of ZTD among the ACs with regard to an unweighted mean is 4.6?mm. The ratio of the analysis noise to the observation noise assessed by the operator/software impact (OSI) model is about 2.5. These and other effects have to be accounted for to improve the intra-technique combination of VLBI-derived troposphere parameters. While the largest systematics caused by inconsistent usage of meteorological data can be avoided and the application of different mapping functions can be considered by applying empirical corrections, the noise has to be modeled in the stochastic model of intra-technique combination. The application of different stochastic models shows no significant effects on the combined parameters but results in different mean formal errors: the mean formal errors of the combined ZTD are 2.3?mm (unweighted), 4.4?mm (diagonal), 8.6?mm [variance component (VC) estimation], and 8.6?mm (operator/software impact, OSI). On the one hand, the OSI model, i.e. the inclusion of off-diagonal elements in the cofactor-matrix, considers the reapplication of observations yielding a factor of about two for mean formal errors as compared to the diagonal approach. On the other hand, the combination based on VC estimation shows large differences among the VCs and exhibits a comparable scaling of formal errors. Thus, for the combination of troposphere parameters a combination of the two extensions of the stochastic model is recommended.     

Synergistic Application of Optical and Radar Data for Archaeological Exploration in the Talakadu Region,Karnataka

Talakadu is a well known historic place situated on bank of the river Cauvery in Mysore district of Karnataka. The place is close to concave side of a prominent meander where large amount of sand has accumulated. It is believed that after construction of a reservoir upstream, sand was exposed to wind action burying the structures of Ganga dynasty and other later kingdoms. A number of buried sites have been identified by archaeological excavations conducted so far. Presently the area forms sand dunes with thick plantation cover. Analysis of RADAR data (fine beam RADARSAT and ENVISAT ASAR) led to identifying a hitherto unknown buried channel through the Old Talakadu town adjoining the excavated archaeological sites. The study suggests that RADAR penetration through the plantation canopy seems to have occurred as observed by comparing with corresponding optical data of LISS-IV. Below the canopy, sand and shrubs on top of the channel (topographically low area) are acting as smooth surface providing dark tone on radar imagery. During field validation GPS was extensively used to navigate through the forest canopy and locate the buried channel, excavated archaeological sites as well as other anomalous patterns. Synergistic application of optical (RESOURCESAT-1 LISS-IV and CARTOSAT-1 & 2) and radar (fine beam RADARSAT and ENVISAT ASAR) data led to identifying remote sensing based guides for archaeological exploration. Integration of known archaeological sites with the identified anomalous patterns was done in GIS environment. This study adds on to the knowledge base of the site and compliments already known information and suggested new areas for further archaeological exploration.     

Spatio-Temporal Coherence Based Technique for Near-Real Time Sea-Ice Identification from Scatterometer Data

The identification of sea-ice has frequently been cited as one of the most important tasks for deriving the sea-ice parameters and to avoid erroneous retrieval of wind vector over sea-ice infested oceans using space-borne scatterometer data. Discrimination between sea-ice and ocean is ambiguous under the high wind and/or thin/scattered ice conditions. The pre-launch technique developed for Oceansat-2, utilizes the dual-polarized QuikSCAT scatterometer data by using the spatio-temporal coherence properties of sea ice in addition to backscatter coefficient and the Active Polarization Ratio. Results were compared with the operational sea-ice products from National Snow and Ice Data Center. The threshold API value of −0.025 was found optimum for sea-ice and ocean discrimination. The overall sea-ice identification accuracy achieved was of the order of 95 per cent, ranging from 92.5% (during December in Southern Hemisphere) to 98% (during March in Northern Hemisphere). The applicability of the algorithm for both the Arctic as well as Antarctic makes it suitable for its operational use with the Oceansat-2 scatterometer data.     

Constructing a Near Real‐time Space‐time Cube to Depict Urban Ambient Air Pollution Scenario

This study adopts a near real‐time space‐time cube approach to portray a dynamic urban air pollution scenario across space and time. Originating from time geography, space‐time cubes provide an approach to integrate spatial and temporal air pollution information into a 3D space. The base of the cube represents the variation of air pollution in a 2D geographical space while the height represents time. This way, the changes of pollution over time can be described by the different component layers of the cube from the base up. The diurnal ambient ozone (O₃) pollution in Houston, Texas is modeled in this study using the space‐time air pollution cube. Two methods, land use regression (LUR) modeling and spatial interpolation, were applied to build the hourly component layers for the air pollution cube. It was found that the LUR modeling performed better than the spatial interpolation in predicting air pollution level. With the availability of real‐time air pollution data, this approach can be extended to produce real‐time air pollution cube is for more accurate air pollution measurement across space and time, which can provide important support to studies in epidemiology, health geography, and environmental regulation.     

Surface Ice Flow Velocity and Tide Retrieval of the Amery Ice Shelf using Precise Point Positioning

Five days of continuous GPS observation data were collected in the frontal zone of the Amery ice shelf and subsequently post-processed using precise point position (PPP) technology based on precise orbit and clock products from the International GNSS service. The surface ice flow velocity of the observed point was derived from PPP to be 2.25 m/day toward the northeast with an azimuth of 41°. Major semi-diurnal and diurnal oceanic tide constituents could be recovered from the 5 days of PPP-derived height variations and compared well with a hydrodynamic ocean tide model. The PPP technique can replace double-difference GPS positioning in remote or hostile environments, and be used to retrieve the surface ice flow velocity without any reference station. Furthermore, the solution can be derived epoch-by-epoch with accuracy in the centimeters to decimeter range.     

Reviewers of the Journal of Geodesy for Volume 79

Acknowledgement for Referees

Reviewers of the Journal of Geodesy for Volume 79     

Effect of signal contamination in matched-filter detection of the signal on a cluttered background

To derive a matched filter for detecting a weak target signal in a hyperspectral image, an estimate of the band-to-band covariance of the target-free background scene is required. We investigate the effects of including some of the target signal in the background scene. Although the covariance is contaminated by the presence of a target signal (there is increased variance in the direction of the target signature), we find that the matched filter is not necessarily affected. In fact, if the variation in plume strength is strictly uncorrelated with the variation in background spectra, the matched filter and its signal-to-clutter ratio (SCR) performance will not be impaired. While there is little a priori reason to expect significant correlation between the plume and the background, there usually is some residual correlation, and this correlation leads to a suppressing effect that limits the SCR obtainable even for strong plumes. These effects are described and quantified analytically, and the crucial role of this correlation is illustrated with some numerical examples using simulated plumes superimposed on real hyperspectral imagery. In one example, we observe an order-of-magnitude loss in SCR for a matched filter based on the contaminated covariance.