A Meissl-modified Vaníček and Kleusberg kernel to reduce the truncation error in gravimetric geoid computations

A deterministic modification of Stokes's integration kernel is presented which reduces the truncation error when regional gravity data are used in conjunction with a global geopotential model to compute a gravimetric geoid. The modification makes use of a combination of two existing modifications from Vaníček and Kleusberg and Meissl. The former modification applies a root mean square minimisation to the upper bound of the truncation error, whilst the latter causes the Fourier series expansion of the truncation error to coverage to zero more rapidly by setting the kernel to zero at the truncation radius. Green's second identity is used to demonstrate that the truncation error converges to zero faster when a Meissl-type modification is made to the Vaníček and Kleusberg kernel. A special case of this modification is proposed by choosing the degree of modification and integration cap-size such that the Vaníček and Kleusberg kernel passes through zero at the truncation radius. Received: 14 October 1996 / Accepted: 20 October 1997     

Spatial–temporal variations of water vapor content over Ethiopia: a study using GPS observations and the ECMWF model

GPS Solutions - We characterize the spatial–temporal variability of integrated water vapor (IWV) in Ethiopia from a network of global positioning system (GPS) stations and the European Center...     

A sun–crown–sensor model and adapted C-correction logic for topographic correction of high resolution forest imagery

Canopy shadowing mediated by topography is an important source of radiometric distortion on remote sensing images of rugged terrain. Topographic correction based on the sun–canopy–sensor (SCS) model significantly improved over those based on the sun–terrain–sensor (STS) model for surfaces with high forest canopy cover, because the SCS model considers and preserves the geotropic nature of trees. The SCS model accounts for sub-pixel canopy shadowing effects and normalizes the sunlit canopy area within a pixel. However, it does not account for mutual shadowing between neighboring pixels. Pixel-to-pixel shadowing is especially apparent for fine resolution satellite images in which individual tree crowns are resolved. This paper proposes a new topographic correction model: the sun–crown–sensor (SCnS) model based on high-resolution satellite imagery (IKONOS) and high-precision LiDAR digital elevation model. An improvement on the C-correction logic with a radiance partitioning method to address the effects of diffuse irradiance is also introduced (SCnS + C). In addition, we incorporate a weighting variable, based on pixel shadow fraction, on the direct and diffuse radiance portions to enhance the retrieval of at-sensor radiance and reflectance of highly shadowed tree pixels and form another variety of SCnS model (SCnS + W). Model evaluation with IKONOS test data showed that the new SCnS model outperformed the STS and SCS models in quantifying the correlation between terrain-regulated illumination factor and at-sensor radiance. Our adapted C-correction logic based on the sun–crown–sensor geometry and radiance partitioning better represented the general additive effects of diffuse radiation than C parameters derived from the STS or SCS models. The weighting factor W_t also significantly enhanced correction results by reducing within-class standard deviation and balancing the mean pixel radiance between sunlit and shaded slopes. We analyzed these improvements with model comparison on the red and near infrared bands. The advantages of SCnS + C and SCnS + W on both bands are expected to facilitate forest classification and change detection applications.     

An appraisal of Indonesia’s immense peat carbon stock using national peatland maps: uncertainties and potential losses from conversion

Background

A large proportion of the world’s tropical peatlands occur in Indonesia where rapid conversion and associated losses of carbon, biodiversity and ecosystem services have brought peatland management to the forefront of Indonesia’s climate mitigation efforts. We evaluated peat volume from two commonly referenced maps of peat distribution and depth published by Wetlands International (WI) and the Indonesian Ministry of Agriculture (MoA), and used regionally specific values of carbon density to calculate carbon stocks.

Results

Peatland extent and volume published in the MoA maps are lower than those in the WI maps, resulting in lower estimates of carbon storage. We estimate Indonesia’s total peat carbon store to be within 13.6 GtC (the low MoA map estimate) and 40.5 GtC (the high WI map estimate) with a best estimate of 28.1 GtC: the midpoint of medium carbon stock estimates derived from WI (30.8 GtC) and MoA (25.3 GtC) maps. This estimate is about half of previous assessments which used an assumed average value of peat thickness for all Indonesian peatlands, and revises the current global tropical peat carbon pool to 75 GtC. Yet, these results do not diminish the significance of Indonesia’s peatlands, which store an estimated 30% more carbon than the biomass of all Indonesian forests. The largest discrepancy between maps is for the Papua province, which accounts for 62–71% of the overall differences in peat area, volume and carbon storage. According to the MoA map, 80% of Indonesian peatlands are <300 cm thick and thus vulnerable to conversion outside of protected areas according to environmental regulations. The carbon contained in these shallower peatlands is conservatively estimated to be 10.6 GtC, equivalent to 42% of Indonesia’s total peat carbon and about 12 years of global emissions from land use change at current rates.

Conclusions

Considering the high uncertainties in peatland extent, volume and carbon storage revealed in this assessment of current maps, a systematic revision of Indonesia’s peat maps to produce a single geospatial reference that is universally accepted would improve national peat carbon storage estimates and greatly benefit carbon cycle research, land use management and spatial planning.

     

Minimum mean squared error (MSE) adjustment and the optimal Tykhonov–Phillips regularization parameter via reproducing best invariant quadratic uniformly unbiased estimates (repro-BIQUUE)

In a linear Gauss–Markov model, the parameter estimates from BLUUE (Best Linear Uniformly Unbiased Estimate) are not robust against possible outliers in the observations. Moreover, by giving up the unbiasedness constraint, the mean squared error (MSE) risk may be further reduced, in particular when the problem is ill-posed. In this paper, the α-weighted S-homBLE (Best homogeneously Linear Estimate) is derived via formulas originally used for variance component estimation on the basis of the repro-BIQUUE (reproducing Best Invariant Quadratic Uniformly Unbiased Estimate) principle in a model with stochastic prior information. In the present model, however, such prior information is not included, which allows the comparison of the stochastic approach (α-weighted S-homBLE) with the well-established algebraic approach of Tykhonov–Phillips regularization, also known as R-HAPS (Hybrid APproximation Solution), whenever the inverse of the “substitute matrix” S exists and is chosen as the R matrix that defines the relative impact of the regularizing term on the final result. The delay in publishing this paper is due to a number of unfortunate complications. It was first submitted as a multi-author paper in two parts. Due to some miscommunication among the original authors, it was reassigned to one of the J Geod special issues, but later reassigned at this author’s request to a standard issue of J Geod. This compounded with a difficulty to find willing reviewers to slow the process. We apologize to the author.     

Towards long-multitemporal change detection using SVI differencing by integrated DWT–ISOCLUS: a model for forest temporal dynamics mapping

Abstract

Characterisation and mapping of land cover/land use within forest areas over long-multitemporal intervals is a complex task. This complexity is mainly due to the location and extent of such areas and, as a consequence, to the lack of full continuous cloud-free coverage of those large regions by one single remote sensing instrument. In order to provide improved long-multitemporal forest change detection using Landsat MSS and ETM + in part of Mt. Kenya rainforest, and to develop a model for forest change monitoring, wavelet transforms analysis was tested against the ISOCLUS algorithm for the derivation of changes in natural forest cover, as determined using four simple ratio-based Vegetation Indices: Simple Ratio (SR), Normalised Difference Vegetation Index (NDVI), Renormalised Difference Vegetation Index (RDVI) and modified simple ratio (MSR). Based on statistical and empirical accuracy assessments, RDVI presented the optimal index for the case study. The overall accuracy statistic of the wavelet derived change/no-change was used to rank the performances of the indices as: RDVI (91.68%), MSR (82.55%), NDVI (79.73%) and SR (65.34%). The integrated discrete wavelet transform–ISOCLUS (DWT–ISOCLUS) result was 42.65% higher than the independent ISOCLUS approach in mapping the change/no-change information. The methodology suggested in this study presents a cost-effective and practical method to detect land-cover changes in support of decision-making for updating forest databases, and for long-term monitoring of vegetation changes from multisensor imagery. The current research contributes to Digital Earth with regards to geo-data acquisition, data mining and representation of one forest systems.     

The welfare benefit of a home’s location: an empirical comparison of spatial and non-spatial model estimates

This article compares results from non-spatial and new spatial methods to examine the reliability of welfare estimates (direct and multiplier effects) for locational housing attributes in Seattle, WA. In particular, we assess if OLS with spatial fixed effects is able to account for the spatial structure in a way that represents a viable alternative to spatial econometric methods. We find that while OLS with spatial fixed effects accounts for more of the spatial structure than simple OLS, it does not account for all of the spatial structure. It thus does not present a viable alternative to the spatial methods. Similar to existing comparisons between results from non-spatial and established spatial methods, we also find that OLS generates higher coefficient and direct effect estimates for both structural and locational housing characteristics than spatial methods do. OLS with spatial fixed effects is closer to the spatial estimates than OLS without fixed effects but remains higher. Finally, a comparison of the direct effects with locally weighted regression results highlights spatial threshold effects that are missed in the global models. Differences between spatial estimators are almost negligible in this study.     

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.     

Estimating rate uncertainty with maximum likelihood: differences between power-law and flicker–random-walk models

Recent studies have documented that global positioning system (GPS) time series of position estimates have temporal correlations which have been modeled as a combination of power-law and white noise processes. When estimating quantities such as a constant rate from GPS time series data, the estimated uncertainties on these quantities are more realistic when using a noise model that includes temporal correlations than simply assuming temporally uncorrelated noise. However, the choice of the specific representation of correlated noise can affect the estimate of uncertainty. For many GPS time series, the background noise can be represented by either: (1) a sum of flicker and random-walk noise or, (2) as a power-law noise model that represents an average of the flicker and random-walk noise. For instance, if the underlying noise model is a combination of flicker and random-walk noise, then incorrectly choosing the power-law model could underestimate the rate uncertainty by a factor of two. Distinguishing between the two alternate noise models is difficult since the flicker component can dominate the assessment of the noise properties because it is spread over a significant portion of the measurable frequency band. But, although not necessarily detectable, the random-walk component can be a major constituent of the estimated rate uncertainty. None the less, it is possible to determine the upper bound on the random-walk noise.     

Contemporary deformation in the Kashmir–Himachal,Garhwal and Kumaon Himalaya: significant insights from 1995–2008 GPS time series

We present new insights on the time-averaged surface velocities, convergence and extension rates along arc-normal transects in Kumaon, Garhwal and Kashmir–Himachal regions in the Indian Himalaya from 13 years of high-precision Global Positioning System (GPS) time series (1995–2008) derived from GPS data at 14 GPS permanent and 42 campaign stations between $29.5{-}35^{\circ }\hbox {N}$ and $76{-}81^{\circ }\hbox {E}$ . The GPS surface horizontal velocities vary significantly from the Higher to Lesser Himalaya and are of the order of 30 to 48 mm/year NE in ITRF 2005 reference frame, and 17 to 2 mm/year SW in an India fixed reference frame indicating that this region is accommodating less than 2 cm/year of the India–Eurasia plate motion ( ${\sim }4~\hbox {cm/year}$ ). The total arc-normal shortening varies between ${\sim }10{-}14~\hbox {mm/year}$ along the different transects of the northwest Himalayan wedge, between the Indo-Tsangpo suture to the north and the Indo-Gangetic foreland to the south indicating high strain accumulation in the Himalayan wedge. This convergence is being accommodated differentially along the arc-normal transects; ${\sim } 5{-}10~\hbox {mm/year}$ in Lesser Himalaya and 3–4 mm/year in Higher Himalaya south of South Tibetan Detachment. Most of the convergence in the Lesser Himalaya of Garhwal and Kumaon is being accommodated just south of the Main Central Thrust fault trace, indicating high strain accumulation in this region which is also consistent with the high seismic activity in this region. In addition, for the first time an arc-normal extension of ${\sim }6~\hbox {mm/year}$ has also been observed in the Tethyan Himalaya of Kumaon. Inverse modeling of GPS-derived surface deformation rates in Garhwal and Kumaon Himalaya using a single dislocation indicate that the Main Himalayan Thrust is locked from the surface to a depth of ${\sim }15{-}20~\hbox {km}$ over a width of 110 km with associated slip rate of ${\sim }16{-}18~\hbox {mm/year}$ . These results indicate that the arc-normal rates in the Northwest Himalaya have a complex deformation pattern involving both convergence and extension, and rigorous seismo-tectonic models in the Himalaya are necessary to account for this pattern. In addition, the results also gave an estimate of co-seismic and post-seismic motion associated with the 1999 Chamoli earthquake, which is modeled to derive the slip and geometry of the rupture plane.     

The static gravity field model DGM-1S from GRACE and GOCE data: computation, validation and an analysis of GOCE mission’s added value

We present a global static model of the Earth’s gravity field entitled DGM-1S based on GRACE and GOCE data. The collection of used data sets includes nearly 7 years of GRACE KBR data and 10 months of GOCE gravity gradient data. The KBR data are transformed with a 3-point differentiation into quantities that are approximately inter-satellite accelerations. Gravity gradients are processed in the instrumental frame. Noise is handled with a frequency-dependent data weighting. DGM-1S is complete to spherical harmonic degree 250 with a Kaula regularization being applied above degree 179. Its performance is compared with a number of other satellite-only GRACE/GOCE models by confronting them with (i) an independent model of the oceanic mean dynamic topography, and (ii) independent KBR and gravity gradient data. The tests reveal a competitive quality for DGM-1S. Importantly, we study added value of GOCE data by comparing the performance of satellite-only GRACE/GOCE models with models produced without GOCE data: either ITG-Grace2010s or EGM2008 depending on which of the two performs better in a given region. The test executed based on independent gravity gradients quantifies this added value as 25–38 % in the continental areas poorly covered with terrestrial gravimetry data (Equatorial Africa, Himalayas, and South America), 7–17 % in those with a good coverage with these data (Australia, North America, and North Eurasia), and 14 % in the oceans. This added value is shown to be almost entirely related to coefficients below degree 200. It is shown that this gain must be entirely attributed to gravity gradients acquired by the mission. The test executed based on an independent model of the mean dynamic topography suggests that problems still seem to exist in satellite-only GRACE/GOCE models over the Pacific ocean, where noticeable deviations between these models and EGM2008 are detected, too.     

Surveying general prospects and challenges of GIS implementation in developing countries: a SWOT–AHP approach

We propose a combined method based on the strengths, weaknesses, opportunities and threats (SWOT) and analytic hierarchy process (AHP) to investigate the challenges and prospects of adopting geographic information systems (GIS) in developing countries. In this context, we identify, group, and analyse SWOT indicators in relation to the main GIS components: data, people, and technology. The relative significance of each SWOT indicator and its related SWOT groups in each GIS component is quantified. The method is then applied in a situation assessment of GIS adoption in the governmental organisations and strategic planning. The SWOT–AHP approach proves to be very useful in identifying and quantifying the relative significance of the major factors affecting GIS implementation, and effectively facilitates GIS strategic planning.      

Coloured Roads on Ordnance Survey First Edition 1 : 2500 Plans and One-Inch Maps 1897–1935, and the Rights Of Way Disclaimer

Abstract

In 1892 a government committee of inquiry into the Ordnance Survey suggested that the Survey should make a distinction between public and private second class roads at the one-inch scale. This study is confined to England and Wales and looks at the nineteenth-century practice of colouring roads on the large-scale plans, and compares the representation of coloured rural roads on early Ordnance Survey one-inch maps with near contemporary highway records of the county of Huntingdonshire, sixteen rural district councils in six counties and an estate in Suffolk. The results show that all but four of the coloured roads on the one-inch maps in these areas are recorded as publicly maintainable highways.     

Cultural and economic residential sorting of Auckland’s population, 1991–2013: an entropy approach

Journal of Geographical Systems - Auckland, the largest city of New Zealand, is one of the most diverse cities in the world, with more than 40% of its population born abroad, more than 200...     

Improvement of the empirical thermospheric model DTM: DTM94 – a comparative review of various temporal variations and prospects in space geodesy applications

The Drag Temperature Model (1978) has been improved to respond better to the actual requirements of space geodesy, especially under extreme solar and geomagnetic conditions. Extended data and an improved algorithm have been considered, leading to valuable improvements. Temporal variations of temperature of the thermopause, total density and major chemical constituent density are reviewed and compared to the DTM94, DTM78 and MSIS86 models. A comparison with data is performed, giving the mean ratios between observed and model values with their root mean squares for different physical and geometric conditions. This comparison is made for the three models and includes data used in the modelling as well as external data. The limits of the thermosphere modelling are discussed. Received: 3 December 1996 / Accepted: 4 August 1997     

Spatial and temporal variations of spatial population accessibility to public hospitals: a case study of rural–urban comparison

Quantification and assessment of nationwide population access to health-care services is a critical undertaking for improving population health and optimizing the performance of national health systems. Rural–urban unbalance of population access to health-care services is widely involved in most of the nations. This unbalance is also potentially affected by varied weather and road conditions. This study investigates the rural and urban performances of public health system by quantifying the spatiotemporal variations of accessibility and assessing the impacts of potential factors. Australian health-care system is used as a case study for the rural–urban comparison of population accessibility. A nationwide travel time-based modified kernel density two-step floating catchment area (MKD2SFCA) model is utilized to compute accessibility of travel time within 30, 60, 120, and 240 min to all public hospitals, hospitals that provide emergency care, and hospitals that provide surgery service, respectively. Results show that accessibility is varied both temporally and spatially, and the rural–urban unbalance is distinct for different types of hospitals. In Australia, from the perspective of spatial distributions of health-care resources, spatial accessibility to all public hospitals in remote and very remote areas is not lower (and may even higher) than that in major cities, but the accessibility to hospitals that provide emergency and surgery services is much higher in major cities than other areas. From the angle of temporal variation of accessibility to public hospitals, reduction of traffic speed is 1.00–3.57% due to precipitation and heavy rain, but it leads to 18–23% and 31–50% of reduction of accessibility in hot-spot and cold-spot regions, respectively, and the impact is severe in New South Wales, Queensland, and Northern Territory during wet seasons. Spatiotemporal analysis for the variations of accessibility can provide quantitative and accurate evidence for geographically local and dynamic strategies of allocation decision-making of medical resources and optimizing health-care systems both locally and nationally.     

Tree Crown Detection,Delineation and Counting in UAV Remote Sensed Images: A Neural Network Based Spectral–Spatial Method

UAVs are fast emerging as a remote sensing platform to complement satellite based remote sensing. Agriculture and ecology is one of the important applications of UAV remote sensing, also known as low altitude remote sensing (LARS). This work demonstrates the use and potential of LARS in agriculture, particularly small holder open field agriculture. Two UAVs are used for remote sensing. The first UAV is a fixed wing aircraft with a high spatial resolution visible spectrum also known as RGB camera as a payload. The second UAV is a quadrotor UAV with an RGB camera interfaced to an on-board single board computer as the payload. LARS was carried out to acquire aerial high spatial resolution RGB images of different farms. Spectral–spatial classification of high spatial resolution RGB images for detection, delineation and counting of tree crowns in the image is presented. Supervised classification is carried out using extreme learning machine (ELM), a single hidden layer feed forward network neural network classifier. ELM was modelled for RGB values as input feature vectors and binary (tree and non-tree pixels) output class. Due to similarities in spectral intensities, some of the non-tree pixels were classified as tree pixels and in order to remove them, spatial classification was performed on the image. Spatial classification was carried out using thresholded geometrical property filtering techniques. Threshold values chosen for carrying out spatial classification were analysed to obtain optimal values. Finally in the delineation and counting, the connected tree crowns were segmented using Watershed algorithm performed on the image after marking individual tree crowns using Distance Transform method. Five representative UAV images captured at different altitudes with different crowns of banana plant, mango trees and coconut trees were used to demonstrate the performance of the proposed method. The performance was compared with the traditional KMeans spectral–spatial method of clustering. Results and comparison of performance parameters of KMeans spectral–spatial and ELM spectral–spatial classification methods are presented. Results indicate that ELM performed better than KMeans.     

Land use and land cover changes over a century (1914–2007) in the Neyyar River Basin,Kerala: a remote sensing and GIS approach

Abstract

Land use and land cover change, perhaps the most significant anthropogenic disturbance to the environment, mainly due to rapid urbanization/industrialization and large scale agricultural activities. In this paper, an attempt has been made to appraise land use/land cover changes over a century (1914–2007) in the Neyyar River Basin (L=56 km; Area = 483.4 km²) in southern Kerala – a biodiversity hot spot in Peninsular India. In this study, digital remote sensing data of the Indian Remote Sensing satellite series I-D (LISS III, 2006–2007) on 1:50,000 scale, Survey of India (SOI) toposheet of 1914 (1:63,360) and 1967 (1:50,000) have been utilized to map various land use/land cover changes. Maps of different periods have been registered and resampled to similar geographic coordinates using ERDAS Imagine 9.0. The most notable changes include decreases in areas of paddy cultivation, mixed crops, scrub lands and evergreen forests, and increases in built-up areas, rubber plantations, dense mixed forests, and water bodies. Further, large scale exploitation of flood plain mud and river sand have reached menacing proportions leading to bank caving and cut offs at channel bends. Conservation of land and water resources forms an important aspect of ecosystem management in the basin.