Accuracy Assessment of Digital Elevation Model Generated by SAR Stereoscopic Technique Using COSMO-Skymed Data

Radargrammetry technique using the stereoscopic synthetic aperture radar (SAR) images is used for the generation of a digital elevation model (DEM) of a region requires only the amplitude images. SAR stereoscopic technique is analogous to the stereo-photogrammetric technique where the optical stereoscopic images are used for DEM generation. While the advantages of the SAR images are their indifference to atmospheric transparency and solar illumination conditions, the side-looking geometry of the SAR increases the complexity in the SAR stereo analysis. The availability of high spatial and temporal resolution SAR data in recent years has facilitated generation of high-resolution DEM with greater vertical accuracy using radargrammetric technique. In the present study, attempt has been made to generate the DEM of Dehra Dun region, India, from the COSMO-Skymed X-band SAR data-pair acquired at 8 days interval through the radargrammetry technique. Here, radargrammetric orientation approach has been adopted to generate the DEM and various issues and processing steps with the radargrammetry technique have been discussed. The DEM was validated with ground measured elevation values using a differential global positioning system and the root-mean-square error of the DEM was found as 7.3 m. The DEM was compared with the reference DEM of the study area generated from the Cartosat-1 stereo data with a model accuracy of 4 m.     

Study on Sampling Rule and Simplification of LiDAR Point Cloud Based on Terrain Complexity

It is difficult to obtain digital elevation model (DEM) in the mountainous regions. As an emerging technology, Light Detection and Ranging (LiDAR) is an enabling technology. However, the amount of points obtained by LiDAR is huge. When processing LiDAR point cloud, huge data will lead to a rapid decline in data processing speed, so it is necessary to thin LiDAR point cloud. In this paper, a new terrain sampling rule had been built based on the integrated terrain complexity, and then based on the rule a LiDAR point cloud simplification method, which was referred as to TCthin, had been proposed. The TCthin method was evaluated by experiments in which XUthin and Lasthin were selected as the TCthin’s comparative methods. The TCthin’s simplification degree was estimated by the simplification rate value, and the TCthin’s simplification quality was evaluated by Root Mean Square Deviation. The experimental results show that the TCthin method can thin LiDAR point cloud effectively and improve the simplification quality, and at 5 m, 10 m, 30 m scale levels, the TCthin method has a good applicability in the areas with different terrain complexity. This study has theoretical and practical value in sampling theory, thinning LiDAR point cloud, building high-precision DEM and so on.     

Sea-Level Change Over the China Sea and Its Vicinity Derived from 25-Year T/P Series Altimeter Data

The characteristics of sea-level change in the China Sea and its vicinity are studied by combining TOPEX/Poseidon (T/P), Jason-1, Jason-2, and Jason-3 altimeter data. First, the sea-surface height is computed by using monthly data via collinear adjustment, regional selection, and crossover adjustment. The sea-level anomaly (SLA) from October 1992 to July 2017 is calculated based on the difference that is obtained by the value derived from the inverse distance weighting method to interpolate the CNES_CLS15 model value at a normal point. By analyzing the satellite data at the same time in orbit, three mean bias groups over the China Sea and its vicinity are obtained: the difference between T/P and Jason-1 is ??11.76 cm, the difference between Jason-1 and Jason-2 is 9.6 cm, and the difference between Jason-2 and Jason-3 is 2.42 cm. To establish an SLA series for 25 years in the study area, the SLAs are corrected. Mean rate of sea-level rise of the Bohai Sea, Yellow Sea, East China Sea, and South China Sea of 4.87 mm/a, 2.68 mm/a, 2.88 mm/a, and 4.67 mm/a, respectively, is found by analyzing the series of SLAs.     

Surface Reflectance Climate Data Records (CDRs) is a Reliable Landsat ETM+ Source to Study Chlorophyll Content in Pecan Orchards

We evaluated the relationships among three Landsat Enhanced Thematic Mapper (ETM+) datasets, top-of-atmosphere (TOA) reflectance, surface reflectance climate data records (surface reflectance-_CDR) and atmospherically corrected images using Fast Line-of-Sight atmospheric analysis of Spectral Hypercubes model (surface reflectance-_FLAASH) and their linkto pecan foliar chlorophyll content(chl-cont). Foliar chlorophyll content as determined with a SPAD meter, and remotely-sensed data were collected from two mature pecan orchards (one grown in a sandy loam and the other in clay loam soil) during the experimental period. Enhanced vegetation index derived from remotely sensed data was correlated to chl-cont. At both orchards, TOA reflectance was significantly lower than surface reflectance within the 550–2400 nm wavelength range. Reflectance from atmospherically corrected images (surface reflectance-_CDR and surface reflectance-_FLAASH) was similar in the shortwave infrared (SWIR: 1550–1750 and 2080–2350 nm) and statistically different in the visible (350–700 nm). Enhanced vegetation index derived from surface reflectance-_CDR and surface reflectance-_FLAASH had higher correlation with chl-cont than TOA. Accordingly, surface reflectance is an essential prerequisite for using Landsat ETM+  data and TOA reflectance could lead to miss-/or underestimate chl-cont in pecan orchards. Interestingly, the correlation comparisons (Williams t test) between surface reflectance-_CDR and chl-cont was statistically similar to the correlation between chl-cont and commercial atmospheric correction model. Overall, surface reflectance-_CDR, which is freely available from the earth explorer portal, is a reliable atmospherically corrected Landsat ETM+ image source to study foliar chlorophyll content in pecan orchards.     

Comparative Analyses of the Point Cloud Produced by Using Close-Range Photogrammetry and Terrestrial Laser Scanning for Rock Surface

Point cloud produced by using theoretically and practically different techniques is one of the most preferred data types in various engineering applications and projects. The advanced methods to obtain point cloud data in terrestrial studies are close range photogrammetry (CRP) and terrestrial laser scanning (TLS). In the TLS technique, separated from the CRP in terms of system structure, denser point cloud at certain intervals can be produced. However, point clouds can be produced with the help of photographs taken at appropriate conditions depending on the hardware and software technologies. Adequate quality photographs can be obtained by consumer grade digital cameras, and photogrammetric software widely used nowadays provides the generation of point cloud support. The tendency and the desire for the TLS are higher since it constitutes a new area of research. Moreover, it is believed that TLS takes the place of CRP, reviewed as antiquated. In this study that is conducted on rock surfaces located at Istanbul Technical University Ayazaga Campus, whether point cloud produced by means photographs can be used instead of point cloud obtained by laser scanner device is investigated. Study is worked on covers approximately area of 30 m?×?10 m. In order to compare the methods, 2D and 3D analyses as well as accuracy assessment were conducted. 2D analysis is areal-based whereas 3D analysis is volume-based. Analyses results showed that point clouds in both cases are similar to each other and can be used for similar other studies. Also, because the factors affecting the accuracy of the basic data and derived product for both methods are quite variable, it was concluded that it is not appropriate to make a choice regardless of the object of interest and the working conditions.     

Science-based approach for credible accounting of mitigation in managed forests

Background

The credibility and effectiveness of country climate targets under the Paris Agreement requires that, in all greenhouse gas (GHG) sectors, the accounted mitigation outcomes reflect genuine deviations from the type and magnitude of activities generating emissions in the base year or baseline. This is challenging for the forestry sector, as the future net emissions can change irrespective of actual management activities, because of age-related stand dynamics resulting from past management and natural disturbances. The solution implemented under the Kyoto Protocol (2013–2020) was accounting mitigation as deviation from a projected (forward-looking) “forest reference level”, which considered the age-related dynamics but also allowed including the assumed future implementation of approved policies. This caused controversies, as unverifiable counterfactual scenarios with inflated future harvest could lead to credits where no change in management has actually occurred, or conversely, failing to reflect in the accounts a policy-driven increase in net emissions. Instead, here we describe an approach to set reference levels based on the projected continuation of documented historical forest management practice, i.e. reflecting age-related dynamics but not the future impact of policies. We illustrate a possible method to implement this approach at the level of the European Union (EU) using the Carbon Budget Model.

Results

Using EU country data, we show that forest sinks between 2013 and 2016 were greater than that assumed in the 2013–2020 EU reference level under the Kyoto Protocol, which would lead to credits of 110–120 Mt CO₂/year (capped at 70–80 Mt CO₂/year, equivalent to 1.3% of 1990 EU total emissions). By modelling the continuation of management practice documented historically (2000–2009), we show that these credits are mostly due to the inclusion in the reference levels of policy-assumed harvest increases that never materialized. With our proposed approach, harvest is expected to increase (12% in 2030 at EU-level, relative to 2000–2009), but more slowly than in current forest reference levels, and only because of age-related dynamics, i.e. increased growing stocks in maturing forests.

Conclusions

Our science-based approach, compatible with the EU post-2020 climate legislation, helps to ensure that only genuine deviations from the continuation of historically documented forest management practices are accounted toward climate targets, therefore enhancing the consistency and comparability across GHG sectors. It provides flexibility for countries to increase harvest in future reference levels when justified by age-related dynamics. It offers a policy-neutral solution to the polarized debate on forest accounting (especially on bioenergy) and supports the credibility of forest sector mitigation under the Paris Agreement.

     

Estimating urban above ground biomass with multi-scale LiDAR

Background

Urban trees have long been valued for providing ecosystem services (mitigation of the “heat island” effect, suppression of air pollution, etc.); more recently the potential of urban forests to store significant above ground biomass (AGB) has also be recognised. However, urban areas pose particular challenges when assessing AGB due to plasticity of tree form, high species diversity as well as heterogeneous and complex land cover. Remote sensing, in particular light detection and ranging (LiDAR), provide a unique opportunity to assess urban AGB by directly measuring tree structure. In this study, terrestrial LiDAR measurements were used to derive new allometry for the London Borough of Camden, that incorporates the wide range of tree structures typical of an urban setting. Using a wall-to-wall airborne LiDAR dataset, individual trees were then identified across the Borough with a new individual tree detection (ITD) method. The new allometry was subsequently applied to the identified trees, generating a Borough-wide estimate of AGB.

Results

Camden has an estimated median AGB density of 51.6 Mg ha^–1 where maximum AGB density is found in pockets of woodland; terrestrial LiDAR-derived AGB estimates suggest these areas are comparable to temperate and tropical forest. Multiple linear regression of terrestrial LiDAR-derived maximum height and projected crown area explained 93% of variance in tree volume, highlighting the utility of these metrics to characterise diverse tree structure. Locally derived allometry provided accurate estimates of tree volume whereas a Borough-wide allometry tended to overestimate AGB in woodland areas. The new ITD method successfully identified individual trees; however, AGB was underestimated by ≤?25% when compared to terrestrial LiDAR, owing to the inability of ITD to resolve crown overlap. A Monte Carlo uncertainty analysis identified assigning wood density values as the largest source of uncertainty when estimating AGB.

Conclusion

Over the coming century global populations are predicted to become increasingly urbanised, leading to an unprecedented expansion of urban land cover. Urban areas will become more important as carbon sinks and effective tools to assess carbon densities in these areas are therefore required. Using multi-scale LiDAR presents an opportunity to achieve this, providing a spatially explicit map of urban forest structure and AGB.

     

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.