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.     

Revisit the calibration errors on experimental slant total electron content (TEC) determined with GPS

The calibration errors on experimental slant total electron content (TEC) determined with global positioning system (GPS) observations is revisited. Instead of the analysis of the calibration errors on the carrier phase leveled to code ionospheric observable, we focus on the accuracy analysis of the undifferenced ambiguity-fixed carrier phase ionospheric observable determined from a global distribution of permanent receivers. The results achieved are: (1) using data from an entire month within the last solar cycle maximum, the undifferenced ambiguity-fixed carrier phase ionospheric observable is found to be over one order of magnitude more accurate than the carrier phase leveled to code ionospheric observable and the raw code ionospheric observable. The observation error of the undifferenced ambiguity-fixed carrier phase ionospheric observable ranges from 0.05 to 0.11 total electron content unit (TECU) while that of the carrier phase leveled to code and the raw code ionospheric observable is from 0.65 to 1.65 and 3.14 to 7.48 TECU, respectively. (2) The time-varying receiver differential code bias (DCB), which presents clear day boundary discontinuity and intra-day variability pattern, contributes the most part of the observation error. This contribution is assessed by the short-term stability of the between-receiver DCB, which ranges from 0.06 to 0.17 TECU in a single day. (3) The remaining part of the observation errors presents a sidereal time cycle pattern, indicating the effects of the multipath. Further, the magnitude of the remaining part implies that the code multipath effects are much reduced. (4) The intra-day variation of the between-receiver DCB of the collocated stations suggests that estimating DCBs as a daily constant can have a mis-modeling error of at least several tenths of 1 TECU.     

Geographic information science and the analysis of place and health

The representation of place is a key theoretical advancement that geographic information science can offer to improve the understanding of environmental determinants of health, but developing robust computational representations of place requires a substantial departure from conventional notions of geographic representation in GIS. Unlike conventional GIS representations based on either objects or locations, we suggest that place representation should incorporate dynamic subjective, experiential, and relational aspects of place, as the influence of place on health behavior concerns not only the features that can be objectively observed at a particular location but also the environmental perceptions of the individual, as molded by biological, social, and experiential characteristics. In addition, assessments of environmental exposures on health outcomes should focus on individuals' time–activity patterns and microenvironment profiles, which form a potentially unique personalized exposure environment for each individual. Addressing these representational challenges via collaborative research has the potential to advance both geographic information science and health research.     

Numerical solution to the oblique derivative boundary value problem on non-uniform grids above the Earth topography

The paper presents a numerical solution of the oblique derivative boundary value problem on and above the Earth’s topography using the finite volume method (FVM). It introduces a novel method for constructing non-uniform hexahedron 3D grids above the Earth’s surface. It is based on an evolution of a surface, which approximates the Earth’s topography, by mean curvature. To obtain optimal shapes of non-uniform 3D grid, the proposed evolution is accompanied by a tangential redistribution of grid nodes. Afterwards, the Laplace equation is discretized using FVM developed for such a non-uniform grid. The oblique derivative boundary condition is treated as a stationary advection equation, and we derive a new upwind type discretization suitable for non-uniform 3D grids. The discretization of the Laplace equation together with the discretization of the oblique derivative boundary condition leads to a linear system of equations. The solution of this system gives the disturbing potential in the whole computational domain including the Earth’s surface. Numerical experiments aim to show properties and demonstrate efficiency of the developed FVM approach. The first experiments study an experimental order of convergence of the method. Then, a reconstruction of the harmonic function on the Earth’s topography, which is generated from the EGM2008 or EIGEN-6C4 global geopotential model, is presented. The obtained FVM solutions show that refining of the computational grid leads to more precise results. The last experiment deals with local gravity field modelling in Slovakia using terrestrial gravity data. The GNSS-levelling test shows accuracy of the obtained local quasigeoid model.     

Integral transformations of gradiometric data onto a GRACE type of observable

Integral transformations of gravitational gradients onto a Gravity Recovery And Climate Experiment (GRACE) type of observable are derived in this article. The gravitational gradients represent components of the gravitational tensor in the local north-oriented frame. The GRACE type of observable corresponds to a difference between two gravitational vectors as projected onto the line of sight between the two GRACE satellites. In total, three integral transformations relating vertical–vertical, vertical–horizontal and horizontal–horizontal gravitational gradients with the GRACE type of observable are provided. Spectral and closed forms of corresponding isotropic kernels are derived for each transformation. Special cases show that the integral transformations are general and relate gravitational gradients to many other quantities of the gravitational field, such as the gravitational vector, and its radial and tangential components. Correctness of the mathematical derivations is validated in a closed-loop simulation using synthetic data.     

Similarity measurement on human mobility data with spatially weighted structural similarity index (SpSSIM)

Understanding diverse characteristics of human mobility provides profound knowledge of urban dynamics and complexity. Human movements are recorded in a variety of data sources and each describes unique mobility characteristics. Revealing similarity and difference in mobility data sources facilitates grasping comprehensive human mobility patterns. This study introduces a new method to measure similarities on two origin–destination (OD) matrices by spatially extending an image‐assessment tool, the structural similarity index (SSIM). The new measurement, spatially weighted SSIM (SpSSIM), utilizes weight matrices to overcome the SSIM sensitivity issue due to the ordering of OD pairs by explicitly defining spatial adjacency. To evaluate SpSSIM, we compared performances between SSIM and SpSSIM with resampling the orders of OD pairs and conducted bootstrapping to test the statistical significance of SpSSIM. As a case study, we compared OD matrices generated from three data sources in San Diego County, CA: U.S. Census‐based Longitudinal Employer–Household Dynamics Origin–Destination employment statistics, Twitter, and Instagram. The case study demonstrated that SpSSIM was able to capture similarities of mobility patterns between datasets that varied by distance. Some regions showed local dissimilarity while the global index indicated they were similar. The results enhance the understanding of complex mobility patterns from various datasets, including social media.     

A hyperspectral experiment over tropical forests based on the EO-1 orbit change and PROSAIL simulation

ABSTRACT

We designed a unique hyperspectral experiment from the Earth Observing One (EO-1) orbit change to evaluate solar illumination effects over tropical forests in Brazil. Ten nadir-viewing Hyperion images collected over a fixed site and period of the year (July to August) were selected for analysis. We evaluated variations in reflectance and in 16 narrowband vegetation indices (VIs) with increasing solar zenith angle (SZA) from the pre-drift (2004–2008) to the EO-1 drift period (2011–2016). To detect changes in reflectance and shadows, we applied spectral mixture analysis (SMA) and principal component analysis (PCA) and calculated the similarity spectral angle (θ) between the vegetation spectra measured with variable SZA. The magnitude of the illumination effects was also evaluated from change-point analysis and nonparametric Mann-Whitney U tests applied over the time series. Finally, we complemented our experiment using the PROSAIL model to simulate the VIs variation with increasing SZA resultant from satellite drift. The results showed significant changes in Hyperion reflectance and VIs, especially when the EO-1 crossed the study area at earlier times and larger SZA in 2015 (9:05 a.m.; SZA = 59°) and 2016 (8:30 a.m.; SZA = 67°). Compared to the pre-drift period (10:30 a.m.; SZA = 45°), the SZA differences of 14° (2015) and 22° (2016) increased the shade fractions and decreased the vegetation brightness. PCA separated the pre-drift and drift reflectance datasets, showing shifts in scores due to changes in brightness. θ increased with SZA, indicating changes in the shape of the vegetation spectra with drift. For most VIs, the change-point analysis indicated 2015 (SZA = 59°) as the predominant year of detected changes. Compared to the EO-1 original orbit, the Plant Senescence Reflectance Index (PSRI), Anthocyanin Reflectance Index (ARI) and Structure Insensitive Pigment Index (SIPI) presented the largest positive changes during drift, while the Photochemical Reflectance Index (PRI), Visible Atmospherically Resistant Index (VARI) and Enhanced Vegetation Index (EVI) had the largest negative changes. The effect size of the illumination geometry on these VIs was large, as indicated by increasing values of the Cohen’s r metric toward 2016. The anisotropy of the Hyperion VIs was generally consistent with that from PROSAIL in the simulated pre-drift and drift periods. Focusing on structural indices, it affected the relationships between VIs and simulated leaf area index (LAI) at large SZA.