From data to decision-making – experiences in Integrated Earth Observation (IEO)

Within the confines of Earth Observation Scientific Knowledge and Technology Transfer in Hungary (EKAT) titled ESA PECS project we gained forum for utilization of our earth observation (EO) knowledge and possibility for further development of experiences. The project aims the horizontal preparation of Hungary for ESA membership – for a dynamically developing sphere in the common part of science, technology and business. For the efficient transfer and dissemination, we preconceived target groups and principles to reach main take-up points and to open new forums of information exchange. Nowadays, user-friendly service, complex supply and real-validation process all have great importance. Joining units of suitable data choice, entire ingestion–processing chains and ergonomic interfaces can only offer a timely and professional solution for the challenges of industry and business sectors. It has to be clearly seen, that information and communication technologies (ICT) and information control becomes a more and more significant part of EO services.     

Monitoring the regional and temporal variability of winter snowfall in the arid Andes using digital NOAA/AVHRR data

Abstract

This paper deals with the spatial distribution and the temporal variability of snowfall in the most arid part of the Andes (18°‐ 28°S) during southern hemisphere winter (May‐September). As the official precipitation data is of poor quality, analyses were carried out by means of digital image processing techniques, using NOAA/AVHRR satellite‐data. Through analysis of 24 different snowfall events from six winters, a previously unknown spatial and temporal precipitation pattern in this remote and unexplored area was revealed. Snowfall is most abundant in the southernmost part of the research area and on the western side of the Andes, indicating the Pacific origin of the snowfall.

Nevertheless, the typical snowfall pattern is modified during different periods of the winter. Three typical time periods could be defined and distinguished from one another. Each of these three periods is characterized by typical weather conditions (cold fronts and “cut‐offs “) leading to a distinct snowfall pattern.

As this study is part of a broader paleoclimatic project, the results will serve as a basis for paleoclimatic reconstruction of past climate. Only by knowing the modern circulation and precipitation patterns is it possible to interpret paleoclimatic signals and archives found in the study area (e.g. paleosol, moraines) correctly.     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Ocean tide loading (OTL) displacements from global and local grids: comparisons to GPS estimates over the shelf of Brittany, France

In this paper we examine OTL displacements detected by GPS stations of a dedicated campaign and validate ocean tide models. Our area of study is the continental shelf of Brittany and Cotentin in France. Brittany is one of the few places in the world where tides provoke loading displacements of ∼10–12 cm vertically and a few cm horizontally. Ocean tide models suffer from important discrepancies in this region. Seven global and regional ocean tide models were tested: FES2004 corrected for K2, TPXO.7.0, TPXO.6.2, GOT00.2, CSR4.0, NAO.99b and the most recent regional grids of the North East Atlantic (NEA2004). These gridded amplitudes and phases of ocean tides were convolved in order to get the predicted OTL displacements using two different algorithms. Data over a period of 3.5 months of 8 GPS campaign stations located on the north coast of Brittany are used, in order to evaluate the geographical distribution of the OTL effect. We have modified and implemented new algorithms in our GPS software, GINS 7.1. GPS OTL constituents are estimated based on 1-day batch solutions. We compare the observed GPS OTL constituents of M₂, S₂, N₂ and K₁ waves with the selected ocean tide models on global and regional grids. Large phase-lag and amplitude discrepancies over 20° and 1.5 cm in the vertical direction in the semi-diurnal band of M₂ between predictions and GPS/models are detected in the Bay of Mont St-Michel. From a least squares spectral analysis of the GPS time-series, significant harmonic peaks in the integer multiples of the orbital periods of the GPS satellites are observed, indicating the existence of multipath effects in the GPS OTL constituents. The GPS OTL observations agree best with FES2004, NEA2004, GOT00.2 and CSR4.0 tide models.     

Contextual Spatiospectral Postreconstruction of Cloud-Contaminated Images

A general method has been proposed recently for the contextual reconstruction of cloud-contaminated areas in multitemporal multispectral images. It is based on the idea of making the prediction process learn from information available in the cloud-free neighborhood of contaminated areas. Though promising, this method does not fully exploit all available information, thus leaving room for further methodological enhancements. This letter presents a postreconstruction methodology for improving the contextual reconstruction process by opportunely capturing spatial and spectral correlations characterizing the considered image. In addition, we propose a solution to a problem that has not yet been addressed in the remote sensing literature, i.e., the generation of an error map beside the reconstructed images to provide end-users with helpful indications about reconstruction reliability. Thorough experiments conducted on a multitemporal sequence of Landsat-7 ETM+ images are reported and discussed.     

GPS observations of the ionospheric F2-layer behavior during the 20th November 2003 geomagnetic storm over South Korea

The ionospheric F2-layer peak density (NmF2) and its height (hmF2) are of great influence on the shape of the ionospheric electron density profile Ne (h) and may be indicative of other physical processes within the ionosphere, especially those due to geomagnetic storms. Such parameters are often estimated using models such as the semiempirical international reference ionosphere (IRI) models or are measured using moderately priced to expensive instrumentation, such as ionosondes or incoherent scatter radars. Global positioning system (GPS) observations have become a powerful tool for mapping high-resolution ionospheric structures, which can be used to study the ionospheric response to geomagnetic storms. In this paper, we describe how 3-D ionospheric electron density profiles were produced from data of the dense permanent Korean GPS network using the tomography reconstruction technique. These profiles are verified by independent ionosonde data. The responses of GPS-derived parameters at the ionospheric F2-layer to the 20th November 2003 geomagnetic storm over South Korea are investigated. A fairly large increase in the electron density at the F2-layer peak (the NmF2) (positive storm) has been observed during this storm, which is accompanied by a significant uplift in the height of the F2 layer peak (the hmF2). This is confirmed by independent ionosonde observations. We suggest that the F2-layer peak height uplift and NmF2 increase are mainly associated with a strong eastward electric field, and are not associated with the increase of the O/N₂ ratio obtained from the GUVI instruments aboard the TIMED satellite. It is also inferred that the increase in NmF2 is not caused by the changes in neutral composition, but is related to other nonchemical effects, such as dynamical changes of vertical ion motions induced by winds and E × B drifts, tides and waves in the mesosphere/lower thermosphere region, which can be dynamically coupled upward to generate ionospheric perturbations and oscillations.     

Seasonal Variation of MODIS Vegetation Indexes and Their Statistical Relationship With Climate Over the Subtropic Evergreen Forest in Zhejiang, China

Moderate Resolution Imaging Spectroradiometer (MODIS) 16-day 1-km vegetation index products, daily temperature, photosynthetically active radiation (PAR), and precipitation from 2001 to 2004 were utilized to analyze the temporal variations of the MODIS normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI), as well as their correlations with climate over the evergreen forested sites in Zhejiang-a humid subtropical region in the southeast of China. The results showed that both NDVI and EVI could discern the seasonal variation of the evergreen forests. Attributed to the sufficient precipitation in the study area, the growth of vegetation is mainly controlled by energy; as a result, NDVI, and especially EVI, is more correlated with temperature and PAR than precipitation. Compared with NDVI, EVI is more sensitive to climate condition and is a better indicator to study vegetation variations in the study region     

A bundle adjustment approach with inner constraints for the scaled orthographic projection

Bundle adjustment is a method for simultaneously calculating both the interior and exterior orientation parameters of a set of images, and the object-space coordinates of the observed points. In the case of long focal length lenses and narrow field-of-view (FOV) imaging situations, collinearity based (perspective projection) algorithms may result in linear dependencies between parameters that cause solution instability. The use of a scaled orthographic projection model based on linear algebraic formulations was therefore adopted to reduce this risk. Using quaternions, a new mathematical model is derived that includes the partial derivatives as well as the inner constraint equations for a scaled orthographic bundle adjustment. The model was then tested using two image sets of a single, small vessel (about 6 m length) with a cube target of known dimensions at two distinct ranges; perspective solutions were also calculated for comparison. RMS residual errors of 0.74-0.78 pixels associated with the new method compare favorably to a residual error range of 0.59-0.74 pixels using a perspective bundle adjustment of the same target points. Relative precisions (as a ratio of target size) of between 1:1650 and 1:750 have been achieved at ranges of 375 m and 850 m, respectively, given comparisons with the known cube dimensions. A third image dataset consisting of a network of 16 images was solved with a 1:2200 relative precision showing the new method can successfully handle high redundancy. For the experiments that were conducted, the new method was found to produce less precise results than the perspective bundle solution for a FOV of 0.50-0.65° where the object fills 5-8% of the image. However, it was found to match the precision of the perspective model (with an uncalibrated camera) for a FOV of 0.20-0.30° where the object of interest fills only 1-2% of the full image.     

Transverse Mercator with an accuracy of a few nanometers

Implementations of two algorithms for the transverse Mercator projection are described; these achieve accuracies close to machine precision. One is based on the exact equations of Thompson and Lee and the other uses an extension of Krüger’s series for the mapping to higher order. The exact method provides an accuracy of 9 nm over the entire ellipsoid, while the errors in the series method are less than 5 nm within 3900 km of the central meridian. In each case, the meridian convergence and scale are also computed with similar accuracy. The speed of the series method is competitive with other less accurate algorithms and the exact method is about five times slower.     

Terrestrial laser scanning to estimate plot-level forest canopy fuel properties

This paper evaluates the potential of a terrestrial laser scanner (TLS) to characterize forest canopy fuel characteristics at plot level. Several canopy properties, namely canopy height, canopy cover, canopy base height and fuel strata gap were estimated. Different approaches were tested to avoid the effect of canopy shadowing on canopy height estimation caused by deployment of the TLS below the canopy. Estimation of canopy height using a grid approach provided a coefficient of determination of R² = 0.81 and an RMSE of 2.47 m. A similar RMSE was obtained using the 99th percentile of the height distribution of the highest points, representing the 1% of the data, although the coefficient of determination was lower (R² = 0.70). Canopy cover (CC) was estimated as a function of the occupied cells of a grid superimposed upon the TLS point clouds. It was found that CC estimates were dependent on the cell size selected, with 3 cm being the optimum resolution for this study. The effect of the zenith view angle on CC estimates was also analyzed. A simple method was developed to estimate canopy base height from the vegetation vertical profiles derived from an occupied/non-occupied voxels approach. Canopy base height was estimated with an RMSE of 3.09 m and an R² = 0.86. Terrestrial laser scanning also provides a unique opportunity to estimate the fuel strata gap (FSG), which has not been previously derived from remotely sensed data. The FSG was also derived from the vegetation vertical profile with an RMSE of 1.53 m and an R² = 0.87.